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1 Introduction 


The discovery of the accelerating expansion of the Universe through observations of distant 
supernovae DEI. together with the theory of inflation |3l IH [5l [6] , suggest the possibility 
that our universe approaches de Sitter geometries in both the far past and the far future. 
However, what is the nature of dark energy and inflation is still one of the deepest puzzles 
in modern physics. We believe that the quantum theory of gravity need to be understood 
before solving this puzzle. 

In an asymptotically flat spacetime a unitary and analytic S-matrix can be well-dehned. 
Even though there is no notion of an S-matrix in asymptotically anti de Sitter spacetime, 
one has the correlation function for the graviton which is equivalent to the correlation 
function of the stress tensor in the conformal held theory on the boundary [7]. Such a 
correspondence between gravity and conformal held theory is crucial for us to understand 
the quantum theory of gravity in the anti de Sitter space. On the other hand, tensor per¬ 
turbations are also generated during inhation [8], whose direct detection would be taken 
as the evidence for inhation, and the non-Gaussianities of gravitational waves is a key 
feature of tensor perturbations as well as the scalar perturbations [H]. The three-point 
correlation function of graviton in de Sitter background has been calculated in [10]. In 
this paper we will investigate the general properties of n-point correlation function of 
graviton, and explicitly compute the four-point function of graviton in de Sitter back¬ 
ground. Even though these gravity correlation functions appear to be outside the reach 
of the experiments occurring in the near future, we hope that the correlation functions of 
graviton in de Sitter background may also shed light on the quantum gravity theory in 
de Sitter space. 

String theory is the most promising quantum theory of gravity. In principle, the 
diagram of the four-point function of graviton at tree level in string theory is unique (see 
Fig.g. Unfortunately, string theory cannot be well formulated in de Sitter space. So we 
do not know how to calculate the four-point function of graviton in de Sitter background 
in string theory. Up to now, general relativity (GR) proposed by Albert Einstein one 
hundred years ago is still the simplest theory of gravity which is consistent with all of 
the experimental data, and every consistent string theory must contain a massless spin-2 
state whose interactions reduce at low energy to general relativity. Here we will expand 
the Einstein-Hilbert action with a positive cosmological constant up to the fourth order 
in terms of gravity fluctuations, and utilize the in-in formalism to calculate the four-point 
correlation function of graviton, including both the contact and exchange diagrams, in de 
Sitter background. 

This paper is organized as follows. In section we expand the Einstein-Hilbert action 
with a positive cosmological constant up to the fourth order in terms of gravitational 
waves. The parity of tensor perturbations and the general properties of n-point correlation 
function of graviton will be discussed in section The four-point function of graviton in 
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Figure 1: The unique diagram of the four-point function of graviton at tree level in string 
theory. 

de Sitter background will be computed explicitly in section Summary and discussion 
will be given in section In addition, a brief review of tensor power spectrum and 
three-point function of graviton in de Sitter background is presented in Appendix and 
we fix the representation of polarization tensors and clarify the general four-momentum 
conhguration into 6 parameters in Appendixj^ The constraints on the correlation function 
of graviton from the de Sitter isometries are discussed in Appendix 


2 The action up to fourth order in terms of tensor 
fluctuations in de Sitter background 

In this paper we focus on the Universe whose dynamics is govern by the Einstein-Hilbert 
action with a positive cosmological constant A, namely 

r 

S=^ dtd^x^iR - 2A), (1) 

where Mpi = is the reduced Planck mass. The background spacetime is de¬ 

scribed by nothing but de Sitter with metric 

= —dt‘^ -I- e^^^dxidxi, (2) 

where H is the Hubble constant related to A by = A/3. 

During a period of exactly de Sitter expansion, there are only two physical degrees 
of freedom corresponding to the tensor fluctuations (gravitational waves). In order to 
make the dynamical degrees of freedom manifest, we call for the the ADM formalism m, 
namely 

= —N‘^dt^ -I- gij{dx^ + N^dt){dx^ + N^dt), (3) 
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( 4 ) 


where N and Ni are lapses and shifts. Now the action Q reads 

r 

S=^ + N-\E,jE^^ - E^)), 

where R^'^\gij\ is the three-dimensional Ricci scalar associated with the metric gtj, and 
the symmetric tensor Eij is dehned by 

~ ~ '^i^j ~ ( 5 ) 

here the dot denotes the derivative with respect to the cosmic time t, and E = g^'^E^j. 
Note that the extrinsic curvature is = N~^Eij. In this sense EijE^^ — E‘^ can be 
taken as the “kinetic term” which governs the dynamics of the extrinsic curvature. In this 
form the spatial coordinate reparameterization is explicitly symmetric, while the time 
reparameterization is not so obvious since iV* mixes up the space and time coordinates. 
In the ADM form the metric gij is dynamical, but N and Ni are Lagrange multipliers. It 
is convenient to choose a gauge for gij in which the time and spatial reparametrizations 
are hxed, and then the metric perturbation can be dehned by 

gij =a^exp{h)-j, ( 6 ) 

where a = and the gravity huctuations are transverse and traceless, i.e. dihij = 0 
and ha = 0. Such a symmetric transverse and traceless metric gij encodes two physical 
degrees of freedom of gravitational waves in GR. 

From Eq. Q, the equations of motion for N and Ni are given by 

/^(3) _ 6^2 _ j^-2 (^EijE^^ - E^) = 0, (7) 

\/i[N-^{E}-ESi)] = 0 . ( 8 ) 

In general, in order to work out the action up to fourth order in terms of gravitational 
waves, one only need to solve the above two equations up to second order in the ADM 
formalism because the third and fourth order perturbations multiplies the hrst and zeroth 
order constraint equation, respectively, which vanishes. See j9] for the discussion in detail. 
Therefore we decompose 


N = 1 + + + ■■■, (9) 

Ni = A^(^) + ATf) + ■ ■ ■. (10) 

Furthermore, the shift can be decomposed into its spin-0 scalar and spin-1 vector 
components in the form of 
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Here a^'^\0{h^). Now one can solve Eqs. M and order by order. Up to 
order of 0{h'^), we get 

„(i) = ^(2) ^ 1^(1) ^ ^(2) ^ ^(1) ^ (^^2) 

+ dkhijdkhij), (13) 

where the prime / denotes the derivative with respect to the conformal time t] = dt'/a{t'). 

Flagging the solntions in (12) and (13) into action (|^ and expanding it np to the 
fonrth order, we obtain 

s = (14) 

where 


^(2) 

^(3) 

^(4) 


Ml, 


Ml, 


dr]d-^xa\h'--h[j - dihijdihij), 


dT]dx^a^{cihikhji - ^hijhki)dkdihij, 


(15) 

(16) 


= Mpi / dr]dx^ 


~a{diij^^^)h'jkdihjk + {h'jkhki - hjkKi) hu 


d^ 


“1“ ^mhij {hjkdrnhkl dm ^jk^kl) ^li ~l~ ^4^ dkhijOihijhkm^lm 
“ 1 ““^^ ^ij ( '^dihjkdyyihkl d-fyihjkdih}^] \ hl^yi 


(17) 


Here the coefficients Cj(= 1) and di{= 1) are nsed to denote the different interaction terms 
in the third and fonrth order action. 


3 Parity of tensor perturbations and the general prop¬ 
erties of n-point correlation function of graviton 


According to the qnadratic action in Eq. (15), the linear pertnrbation equation can be 


solved in the Fourier space. So one can decompose the perturbations into momentum 
modes as follows 

r rj'iU „ 

hij{r],x) = / ——hij(r/,k)e**"'^, (18) 


where 


(27rY 


hijiv.k) = 


(19) 


efj(k) is the polarization tensor which has the following properties 


4(k) = 4(k), eti{k) = kjeUk) = 0. 


( 20 ) 
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For convenience, the polarization tensor is normalized as 

e5(k)e5*(k)=4i„„. (21) 

Usually one can adopt the linear or circular polarization. 


3.1 Circular polarization 


In the literatures the circular polarization is widely used because they are the helicity 
eigenstates of gravitational waves. For the circular polarization, the polarization tensor 
takes the form 


0 


4 ( k ) = 1^-10 
0 0 0 


1 -z 0 

e--^(k) = I -1 0 

0 0 0 


( 22 ) 


where ±2 denote two circularly polarized states which propagate along z direction. From 
the above formula, the circular polarization tensors satisfy 


4(-k) = ey(k) = (4(k))*, 


(23) 


where A = ±2. The most general solution of quantized tensor perturbation can be written 
by 


Kiv) = '0fcaA(k) + V’>l(-k), 
where '0^ is the mode function which satishes 

V’fc - -i^k + = 0, 

V 

and k = |k|. After choosing the standard Bunch-Davies vacuum 
function is given by 

Mv) = • 

Now the creation and annihilation operators are normalized as 


“ii(ki),al,(k2) 


(2T)3i’(k, - kjij.,. 




(24) 

(26) 

the solution of mode 
(26) 

(27) 


From Eq. (24), the operator satishes 


hi = (28) 

On the other hand, since gravity huctuations hij are tensor perturbations, hij follows 

Vhij{r],x.)V~^ = hij{r], -x) . (29) 
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under parity transformation, or equivalently 


Vh^V-^ = hzi, . 


( 30 ) 


Since GR as well as the vacuum state is parity invariant (namely "PlO) = |0)), the u-point 
function of graviton satisfies 


(d; ■■■'"!!;) = • ■ ■ vh^-v-'v) = (vh^v-'^ ■ ■ ■ -p/itP"') 


(31) 


Combining with Eqs. (28) and (30), we hnd 


K\ ■ ■ ■ i4:) = ■ ■ ■ let:) = ('‘d*’ ■ ■ • '‘C) = (Kf ■ ■ ■ kZ)' 


(32) 


which indicates that the n-point function for helicity state possess the symmetry that it 
becomes its complex conjugation when interchanging +2 o —2. 

When all kj lie in the same plane, without loss of generality by assuming that all k,; 
he in y-z plane, there is another symmetry under the mirror reflection about y-z plane. 
The transformation matrix is = diag{—1,1,1}. Denote x = {x,y,z), k = {kx, ky, kz), 
then = {—x,y,z) = x, M^^^k = {—kx,ky,kz) = k, and Myze^{k)My^ = e'^(k). Thus 

under such a mirror reflection. 


or equivalently. 


■^k^yzhij{TI^'xZjJk\y^ {Js/\.yz^ik(^^yz)jlhklij] 1^^ i 


ZZ4(k)MyXtMZ = ZZ4j(k')l4 . 


(33) 

(34) 


Since k is in the y-z plane, = 0, k = k, and then the above equation reads 

ZZeZk)MyJ4M;t = ^e.yk)fti. (35) 

A A 

For any k lies in y-z plane, e'^(k) can be expressed by e'^(k) = Ri(0)e^(z)Ri(d)'^, where 
the definition of Ri(0) is given in (89). Notice that commutes with Ri(0), and hence 

i\k) = M,,e^(k)Mj', = R^(0)Myze\z)M^M0)^ . (36) 


From Fq. (22), one can easily get 


M,,e"(k)M;, = e-'(k), 


(37) 


and then 


e-'(k) = e-^(k) 


(38) 
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Combining with Eq. (35), we find 



( 39 ) 


Therefore, 



(40) 


which indicates that the n-point function for circularly polarized states is always real 
when all momenta are in the same plane. This conclusion is quite nontrivial since is 


not a Hermitian operator for the circularly polarized state, i.e. ^ (^k)^ (28))- 


The bispectrum (three-point function) of graviton is real for the circularly polarized 
states because the three momentum vectors must lie in the same plane due to the energy- 
momentum conservation. But the n(> 4)-point function for the circularly polarized state 
is complex in general. 

3.2 Linear polarization 

The linearly polarized states can be given by a linear combination of circularly polarized 
states, namely 



(41) 


(42) 


and the linear polarization tensors are related to the circular ones by 



(43) 


(44) 


where -|- and x denote two linearly polarized states. One can easily find that the linear 
polarization tensors satisfy 



(45) 


For linear polarization, the most general solution of quantized tensor perturbation be¬ 
comes 


^k(^) = V’fcas(k) + (-l)'2"^/,*al(-k). 


(46) 
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where s = +1, —1 for +, x polarized states respectively. Here ipk is same as that in cir¬ 
cular polarization, and the commutators between the creation and annihilation operators 
are 


Ki(ki).aL(k2)] 


1 


(27r)35='(ki - 


k2)(5sj^s2 • 


(47) 


From Eq. (46), the operator satishes 


K = ('XJ*. K = 


(48) 


Similar to the circular polarization, under the parity transformation the operator h^. 
follows 


Vh*V-^ = Th^V-^ = -h\ . 


(49) 


Therefore 






r 


(50) 


where Eqs. (48) and (49) are taken into account. We conclude that the n-point functions 
of gravitational waves in terms of linear polarized state are always real whether all k* he 
in the same plane or not. Because a real u-point function is convenient in analysis and the 
n-point function in terms of circularly polarized state can be easily derived from n-point 
function in terms of linearly polarized state by using (41), we will compute the four-point 
correlation function of graviton in terms of linearly polarized state in this paper. 


Before closing this subsection, we also want to illustrate some simple properties of 
n-point function of graviton when all k, are in the same plane. Similar to the former 
subsection, assuming that all kj lie in the y-z plane, and under a mirror reflection My^ 
transformation, we have 


M„.e+(k)MT = e+(k). M„.e>'(k)MT = -e^k), 


(51) 


and 


= ftp = -ftp (52) 

Denoting m as the number of x linear polarized modes in the n-point function, we have 

= . (53) 

Thus the n-point function for linearly polarized state equals zero if m is an odd number 
when all the momentum vectors he in the same plan. This property can be seen explicitly 
for the four-point function in the next section. 
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4 The four-point correlation function of graviton in 
de Sitter background 

In this section we will use the in-in formalism to calculate the four-point correlation 
function of graviton during a period of exactly de Sitter expansion, including the contact 
and exchange diagrams. A brief review of tensor power spectrum and the three-point 
function of graviton is given in Appendix 


4.1 In-in formalism 


In this subsection we will introduce the Schwinger-Keldysh in-in formalism [13] to calcu¬ 
late the n-point function of gravity fluctuations during inflation. The expectation value 
of some product Q{t) of held operators at hxed time t is given by 

{Q{t)) = {in\Q{t)\in), (54) 


where \in) is the initial state at very early times when the wavelength is deep inside the 
horizon. In the interaction picture, we have 


m)) = (oi 


*/tn Hl(T')dT' 


Te^o 


Q\t) 


Te 


-*/4 Hl{T'')dT" 


| 0 ) 


(55) 


where Hj is the interaction part of the Hamiltonian. Here T denotes a time-ordered prod¬ 
uct, T denotes an anti-time-ordered product, and is the product Q in the interaction 
picture. Utilizzing the Dyson series, one can get 


iQit)) = 

n=0 




dtn{[Hjitn), [Hj{tn-l), • • • , [Hiih), Q\t)] ■■■]]) 




(56) 


which is called the commutator form. Note that the n-th order in-in formalism have n 
commutators. Supposing that operator Q is Hermitian, i.e. Q = , the hrst and second 

order terms read 

(Q(f))i =2Im [ dh{0\Q\t)Hi{h)\0) , (57) 

Jto 


{Q{t ))2 = I dti I dt 2 { 0 \Hi{ti)Q^{t)Hi{t 2 )\ 0 ) 


Jto 

- 2Re 


'to 


*1 / . 


(58) 


'to 


'to 


Here we need to stress that we should use Eq. (56) to calculate the n-point function if 
the operator Q is not Hermitian. Since there are no time derivative terms in the cubic 
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order interaction Lagrangian for the gravity fluctuations in de Sitter background, the 
interaction part of Hamiltonian follows that Hj = — J dx^Cj up to the fourth order. 

In general the the n-point correlation function of gravitational waves in the de Sitter 
background takes the form 




2=1 


H 


M 


pi 


2(n-l) 


and 


(ki,k2,---,k„) = 


'-J~SiS2--Sn 

(ki,k2, • • • ,k„) , 

(59) 

‘(ki,k2, • 

■ ■) k„) 

5 

(60) 


where iniplies the energy-momentum conservation. Furthermore, the con¬ 

straints on the function of T*^^^'"^”(ki, k 2 , • • •, k„) from the de Sitter isometries are dis¬ 
cussed in Appendix 0 In the next two subsections we will use the in-in formalism to 
explicitly calculate the amplitude ks, k 4 ) which is expected to be different 

for different interaction term. 


4.2 Contact diagram 

In this subsection we will work out the contribution to the four-point correlation function 
of graviton from the contact diagram (see Fig. [^. At the tree level, we only need to use 



Figure 2: Contact diagram. 


the hrst order in-in formalism 

{K\hXK)c 



(61) 
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where = —f d?xC^^\ Remind that Eq. (57) is not applicable in this case because 
the operator is not Hermitian. Before using the above in-in formalism to 

compute four-point function of graviton, we need to hnd the solution of from the 


constraint equation (13). Considering Eq. (18), we can easily hnd 


y^>(x,)() = / (ft' (ki,i;)ft',(k2,i() -ki • k2fty(ki,i;)fty(k2,)() 


(2x)t 


iki2-x 




, (62) 


where kjj = kj -I- kj, kfj = k^ + k'^^ and kij = |kjj|. Now adopting the in-in formalism in 


Eq. (61), after a lengthy but straightforward computation, we get 

X-«(k..k2.k,.k,) = ^.2^.2 


do 


''12 


+ -^Uklkl - klkl)eke%el,ei 


+ X 3 


^(ki ■ kg - ki • k 2 )eke%eliei - 


lm^mn^2 


8" 


+ — I —^k^k^ + /cgfeo" I ehe~i.ei,e 


1 3 4 


23 perms., 
(63) 


where 

Xi = -^{-4{K + 3k4)klk^kl-2ki-k2kl[6kik2k4 + K\2K-k3) 

+2K{kik2 +kiki +k2ki)]}, (64) 

X 2 = [K (K -|- fcg -|- k^) -f 2fc3/c4] , (65) 

X 3 = 2 |iC - ^[2/CiA;2fc3/i^4 + + ^2^3 + ^2^4 + ^ 3 ^ 4 ) 

+K{kik2k3 + kik^ki + kik2k4 + fc2fc3^4))|, (66) 

K = ki + k 2 + ks + ki, and eh = eh(k;). Here the divergent parts from performing the 
time integration are counteracted. In order to sketch out some simple features of the 
four-point function of graviton from the contact diagram, we will consider some special 
conhgurations in the momentum space. 

First of all, we consider the bilateral squeezed limit in which kg = k4 — ~ki2/2, 
ki = k 2 = k and k^ = k^ = k ^ 0. See Fig. In this limit, all k, he in the same plane. 
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Figure 3: Bilateral squeezed limit. 


and there are only four non-zero amplitude, namely 

1 


'y-X XXX _ 

^XX++ _ 

^++xx _ 

J++++ ^ 


96k^k^ 

-=- ( d 2 + 6^3 -|- 86^4) = — - -, 

96k^k^ 96fc3A;6 

1 25 1 

(^2 + 6^3 -|- 36d4 — ISds) = 


25 1 

96 


96k^k^ 

1 


96 fcs^^e’ 


192fc3fc6 


(da - 72^4) = 


71 


192 fc3fc6 


(67) 

( 68 ) 

(69) 

(70) 


We see that the interaction term denoted by di in the fourth order action does not make 
a contribution to the four-point function in the bilateral squeezed limit, and all of the 
non-zero amplitudes have the same scaling behavior (hhhh) oc T oc l/{k^k^) which is 
divergent as 1/fc® in the limit of /c —)■ 0. This is similar to the local-form four-point 
function of scalar perturbation from the contact diagram. 

Secondly, we consider the equilateral quadrangle momentum conhguration in which 
ki = k 2 = k^ = ki = k and = 02 = 0i = 02 = /3. In this conhguration, the two 
angles in Fig. namely a G [0, vr] and 0 G [0,7r/2], are free parameters. Since there 
is not a special momentum mode among these four momentum modes, there are only 
hve independent shapes, i.e. ^ 'y-xxx + ^ -y-xx-i-i-^ -y-x-t-t-n- 7~^+++. The shapes of 

these hve amplitudes show up in Fig. 'y-xxx+ _ -y-x-i-i-i- _ g q, = q or a = tt because 
in these two cases all of four momenta he in the same plane and there are odd number of 
X linearly polarized states. Furthermore, in the limit of 0 —?• 7r/2 corresponding to the 
diamond squeezed limit (see Figj^ in which ki 2 —)■ 0, the four-point function of graviton 
from contact diagram is regular. 


4.3 Exchange diagram 


The third order perturbation Hamiltonian can also make a contribution to the four-point 
function of graviton through the exchange diagram (see Fig. [^. Applying the second 
order in-in formalism, the four-point function from the exchange diagram takes the form 










|o). 


(71) 
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Figure 4: The five independent shapes of trispectrum from contact diagram in the equi¬ 
lateral quadrangle momentum configuration. 
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Figure 6: Exchange diagram, 
where the third order interaction Hamiltonian is given by 

r 1 

^ / dx^a^{cihikhji - -C2hijhki)dkdihij. (72) 

It is remarkable that when doing time double integrals from —cx) to 0 we will get several 
divergent terms due to high negative power of conformal time.Fortunately, these diver¬ 
gences cancel each other exactly when we do permutations of ki, k 2 , ks and k 4 . Finally 
we get the amplitude of the four-point function from the exchange diagram as follows 

^U.,S3S4(k^^k2,k3,k4) (73) 

= “ ^C25*(1,2))(ciH*( 3,4) - ^C25"(3,4)) + 23 perms, 

12 c 
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where 


2^1234 

Bi 

B 2 

Bs 

B, 

B, 

and 


B^ + B2- 2BsB, - 2B,, (74) 

3,2 4), 

-kl, - l(ki + kl + kl + ki) - 2 -h- k,)K - \k^'\ , 

77 + ^ {kik2ki2 + K{ki2{ki + ^ 2 ) - kik2)), 

k^kjku + {ks + k4 + ki2){k3k4 + kujks + ^ 4 )) 

{ki2 + ^3 + kiY 

; A K\ki2 - 7i - k2) + K{ku{ki + ^ 2 ) - 2kik2) + 2ki2kM: 

-^'^(^3 + M + ^ 12 ) 


(A;i 2 + ki + /C 2 ) 


ki2kik2 + {ki2 + /ci + k2){ki2ki + ki2k2 + /C 1 /C 2 ) 
{ki2 + /ci + ^ 2 )^ 


A%i,3) 


+ 7“e:Uk.)A;}e:,(k,,)e:^,(k,) 
+^-e:,(k,,)A;}e:^(k,)e:Kk,) . 


(75) 


Note that i,j = 1, 2, 3,4 represent ith, jth momentum and need not do summation, while 
the indices a, b,c,d = 1 , 2 ,3 are the component indices. 

Analogous to the analysis of contact diagram, in the bilateral squeezed limit, we hnd 

^S1.2S3S4 ^ —c2(si + S2)(l + S3)(l + 54 )^^, (76) 

which indicates that there are only two non-zero components, i.e. and 7 "^^++ 

and only the interaction term denoted by C 2 in the third order action contributes to the 
four-point function in the bilateral squeezed limit. The non-zero components diverge on 
behavior of 1/k^ in the limit of /c —)■ 0 as well. 

Secondly, let’s consider the diamond squeezed limit, namely ki = k 2 = k^ = k/^ = k 
and ki 2 —)■ 0. In this limit, we have 

^SIS2S3S4 ^ Ac2cos2a(si+S 2 )(s 3 + S 4 )p^ • (77) 

Again we see that only the interaction term denoted by C 2 in the third order action makes 
a contribution to the four-point function in this limit. Unlike that from contact diagram, 
the amplitude of four-point function from the exchange diagram goes like I//C 12 which is 
divergent in the limit of ki 2 —t 0. This behavior is similar to the local-form four-point 
function of scalar perturbation from exchange diagram. 
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5 Summary and discussion 


In this paper we dig out some general properties of the n-point function of graviton 
due to the parity, and explicitly compute the four-point correlation function of graviton, 
including both the contact and exchange diagrams, in the de Sitter approximation in 
the framework of GR. One can expect that the the simplicity of the results presented 
in this paper is lost when we go to an inflationary background. However, our results 
are still the leading approximation in the slow-roll expansion. In the framework of GR, 
we hnd {hhhh)/{hh)^ ~ 0{1). Furthermore, the four-point function of graviton in the 
de Sitter background has similar behavior to the local-form four-point function of scalar 
perturbation in some special momentum conhgurations. 

The three-point function of graviton has been extended to the inflationary background 
in [13], generalized G-inflation in [15] and Horava-Lifshitz gravity in [16]. It is worthy 
coming back to the four-point function of graviton in more general scenario in the near 
future. 
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A A brief review of tensor power spectrum and three- 
point correlation function of graviton in de Sitter 
background 


Following the standard technique in quantum held theory, from Eqs. (26), (46) and (47), 
one can easily get 


{fcK'C) = (2>r)"i‘”(k,+k2)i, 


H 

2kl 


(78) 


Gonsidering the normalisation of the polarization tensor and summing over the two degrees 
of freedom of gravitational waves in GR, the amplitude of tensor power spectrum is given 
by Pt = 2{H‘^/MpiY /which is the same as that in the literatures. 

The three-point correlation function of graviton in terms of circularly polarized states 
in de Sitter background has been calculated in [S]. Here we adopt the linearly polarized 
states and distinguish the different contributions from two interaction terms denoted by 


Cl and C 2 in the third order action (16). Using the the hrst order in-in formalism, we 
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obtain 




2=1 




M, 


pi 


where 


^'■««(ki,kj,k3) = - 


\ i<j 


{kik2h] 


/C 1 /C 2 /C 3 


(79) 


/C 2 


X 


{cikYkYiYjYs - ? 44 ^ 3 eL 4 ) + 5 perms. 


(80) 


and 1C = ki + k 2 + k^. The energy-momentum conservation implies that all the momenta 
kj in the three-point function must he in the same plane and then only the amplitudes of 
three-point function with even number of x polarized states do not equal zero. Taking 
the permutations into account, there are only two independent non-zero components 


^xx + (ki,k2,k3) = 


(dci/cg C2{kl + kl- kl)) 
lQ\/2kik2k‘\K, 


—kik2kz — 1C kikj + KC' 


i<j 


X 


^+++(ki,k2,k3) = 


8kik2ks — 4/C kikj -|- /C^ j , 
i<j / 

8kik2kz — 4/C kikj + JC^ 


i<j 


X 


Uy/2{kik2k^flC 

kik2k^ + 1C kikj — 1C 

i<j 

8 c 2 (^ kikjY - 8 ( 2 ci - C 2 )lCkik 2 kz 


i<j 


+4(2ci - 3 c2)/C2 ^ kikj + (3c 2 - 2ci)/C^ 


i<j 


In the squeezed limit, i.e. ki = k 2 = k and k^ —)■ 0, we have 

3c2k^ 


^ 5 iS 253 _ 


16^2 


(Sl + S2)(l + ^ 3 ). 


(81) 


(82) 


(83) 


which says that only the interaction term denoted by C 2 makes contributions to the three- 
point function in the squeezed limit. Note that s = -|-l, —lfor-|-, x polarized states 
respectively. From the above formula, there are only two non-zero components, i.e. 

3c2k^ 


^++- 1 - _ _^XX+ _ 


4^2’ 


(84) 


in the squeezed limit {k^ —)■ 0). It is similar to the local-form three-point function of 
scalar perturbation during inflation. 
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B Polarization tensors 


In this section, we will fix the representation of polarization tensors. Due to the mo¬ 
mentum conservation ki -|- k 2 -|- ks -|- k 4 = 0 for the four-point function, the momentum 
configuration can be described by six parameters, i.e. ki 2 , 9i, 6 * 2 , 0 i, 02 and a, which are 
illustrated in Fig. Here, without loss of generality, we assume that ki and k 2 lie in the 



Figure 7: Momentum configuration. 

x-z plane, and —ki 2 stays along z = ( 0 , 0 , 1 )"*" direction. And then k* can be expressed by 


ki 

= kiR2{n - 9i)z, 

(85) 

k2 

= /c2722(vr-F 0i)z, 

( 86 ) 

ka 

= /c3/23(a)i?2 ( 6 * 2 ) 2 , 

(87) 

k4 

= /C4i?3(7r-f Q!)i?2(02)z, 

( 88 ) 
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where RiiOi) denotes rotating 6i about Xi axis: 


Ri{ei) 


R2{02) 


RsiOs) 


/ 1 0 0 
j 0 cos 6*1 — sin 9i 
y 0 sin 6*1 cos 6*1 

cos 6*2 0 sin 6*2 \ 

0 10, 

— sin 6*2 0 cos 6*2 J 



( cos 6*3 
sin 6*3 
0 


— sin 6*3 
cos 6*3 
0 



Now let’s start with the linear polarization tensor along z 




1+.S 1—s n 

9 9 ^ 

= V2| 


(89) 

(90) 

(91) 


(92) 


where s = 1 (—1) for + (x) polarization state respectively. And then the polarization 
tensors along k* directions are given by 


e«(ki) 

= i?2(7r - 6*i)e^(z)i?2(7r - 6 * 1 )^^, 

(93) 

e*(k2) 

= i?2(7r-6 0i)e®(z)i?2(7r + 0 i)'^. 

(94) 

e*(k3) 

= R3{a)R2{e2)e^{z)R2{e2fR3{af, 

(95) 

e*(k4) 

= Rsiir + a)i?2(02)e''(z)i?2(02)'^7?3(7r -h 

(96) 


C The constraints on the correlation function of gravi¬ 
ton from the de Sitter isometries 

There are 10 isometries in the four-dimensional de Sitter spacetime, including 3 spatial 
rotations, 3 spatial translations, 1 dilatation and 3 special conformal transformations 
(SCT). For de sitter space, the spatial rotations and translations are trivial symmetries 
which just indicates isotropy and homogeneity. 

In the conformal frame, = {—drf -|- dx‘^)/7f and there is a scaling symmetry under 
the transformation rj —)■ Xr], x —)■ Ax. Therefore we have 

(4iii(^i’d)--4:i„(xi,r/)) = (Axi,A?7)...h";jJAx2,Ar/)) , (97) 

which reads 

KM-KHv)} = ( 98 ) 
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for the Fourier modes. From Eq. (59), taking the late time limit (i.e. rj —)■ 0), we hnd 


r- *"(11,, ■■■,!£„) = *"(ki/A,..., k„/A), (99) 

or equivalently 

n 

3(u-l) + ^k„-ak, 

a=l 

The dilatation invariance implies that the spectrum of gravity fluctuations is scale-invariant. 

However, the SCT, namely 77 —)■ r; — 277 (b-x) and x — x-|-b(—77^-1-x^)— 2 x(b-x), where 
b is inhnitesimal, changes the spatial slices and then violates the ADM decomposition. 
Even though in the late time limit SCT becomes 77 —)■ 77 and x — ?■ x -|- b(x^) — 2 x(b ■ x), 
which doesn’t change the spatial slices, the gravity fluctuation hij breaks the traceless 
and transverse conditions under this transformation UTl. In this sense, SCT is not a good 
symmetry for gravitational waves in our gauge-hxing system. Even though SCT may 
provide a constraint on the structure of the n-point function of graviton in the de Sitter 
background, how to work it out is still an open question. 

These four isometries are absent in the realistic inflationary background and the struc¬ 
tures of the 77-point function in the de Sitter background should be corrected up to the 
slow-roll parameters. 


T(ki, ...,k„) = 0. 


( 100 ) 
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